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Abstract. We present a sum-rule extraction of decay constants of heavy mesons from the two-point correlator of heavy-light 
pseudoscalar currents. Our primary concern is to control the uncertainties of the decay constants, induced by both input QCD 
parameters and limited accuracy of the sum-rule method. Gaining this control is possible by applying our novel procedure for 
the extraction of hadron observables utilizing Borel-parameter-depending dual thresholds. For the charmed mesons, we obtain 
fjj = (206.2 ±7. 3 (ope) ±5.1( syst j) MeV and fy s = (245.3 ± 15.7j OPE ) ±4.5( syst )) MeV. In the case of the beauty mesons, the 
decay constants prove to be extremely sensitive to the exact value of the i-quark MS mass mt,(mt,). By matching our sum-rule 
prediction for fg to the lattice outcomes, the very accurate £>-mass value m^{m^) = (4.245 ±0.025) GeV is found, which yields 
f B = (193.4 ± 12.3 ( ope) ±4.3 (syst) ) MeV and fa = (232.5 ± 18.6 (OPE) ±2.4 (syst) ) MeV. 
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QUARK-HADRON DUALITY 

The calculation of the decay constants fp of ground-state 
heavy pseudoscalar mesons P by QCD sum rules [1, 2] is 
a complicated problem: First, a reliable operator product 
expansion (OPE) for the "Borelized" correlation function 
of two pseudoscalar heavy-light currents has to be found. 
Second, even if all the parameters of this OPE are known 
precisely, the knowledge of only a truncated OPE for the 
correlator allows to extract bound-state observables with 
only finite accuracy, reflecting an inherent uncertainty of 
the QCD sum-rule approach. Controlling this uncertainty 
constitutes a delicate problem for actual applications [3]. 

Recall one essential feature of the sum-rule extractions 
of decay constants: the quark-hadron duality assumption 
entails a (merely approximate) relation between hadronic 
ground-state contribution and OPE with the "QCD-level" 
correlator cut at some effective continuum threshold s e ff : 

f^M Q exp(-M Q x) = n dua i(T,s e ff) (1) 

■ v dT 

= J diexp(-ST)p P en(s)+n powei .(<r;). 
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Here, the perturbative spectral density Ppert(i) is obtained 
as a series expansion in powers of the strong coupling a s : 

p peni s)=pM(s) + ^p^(s) + $pW(s) + -. 

Obviously, in order to extract a decay constant /g one has 
to find a way to fix the effective continuum threshold s e ff . 



A crucial albeit very trivial observation is that s e s must 
be a function of r, otherwise the l.h.s. and the r.h.s. of (1) 
would exhibit a different T-behaviour. The exact effective 
continuum threshold — corresponding to the true values 
of hadron mass and decay constant on the l.h.s. of (1) — 
is, of course, not known. Therefore, our idea of extracting 
hadron parameters from sum rules consists in attempting 
(i) to find a reliable approximation to the exact threshold 
5 e ff and (ii) to control the accuracy of this approximation. 
In a recent series of publications [4], we have constructed 
all the associated procedures, techniques and algorithms. 

We define a dual invariant mass Md ua i and a dual decay 
constant /dual (Mq still denoting the true hadron mass) by 

? d 

^dual(T) = -^logn dua i(T,i eff (T)), (2) 

/dlaiO) =M e 4 exp(M^T)n dual (T,s eff (T)). (3) 

In case the ground-state mass Mq is known, the deviation 
of the dual ground-state mass Md ua i from its actual value 
Mq yields an indication of the excited-state contributions 
picked up by our dual correlator. Assuming some specific 
functional shape for our effective threshold and requiring 
least deviation of the dual mass (2) from its actual value 
in the Borel window leads to a variational solution for the 
effective threshold. With i e ff(T) at our disposal we get the 
decay constant from (3). The standard assumption for the 
effective threshold is that it is a (f-independent) constant. 
In addition to such crude approximation we also consider 
polynomials in T. In fact, T-dependent thresholds greatly 
facilitate reproducing the actual mass value. This implies 
that a dual correlator with T-dependent threshold isolates 



the ground state much better and is less contaminated by 
excited states than a dual correlator with the conventional 
T-independent threshold. As consequence, the accuracies 
of extracted hadron observables are drastically improved. 
Recent experience from potential models reveals that the 
band of values obtained from linear, quadratic, and cubic 
Ansatze for the effective threshold encompasses the true 
value of the decay constant [4]. Moreover, we could show 
that the extraction procedures in quantum mechanics and 
in QCD are even quantitatively very similar [5]. Here, we 
report our results [1, 6] for heavy-meson decay constants. 



OPE AND HEAVY-QUARK MASSES 

For heavy-light correlators and emerging decay constants 
the choice of the precise scheme adopted for defining the 
heavy-quark mass has a great impact. We utilize the OPE 
for this correlator to three-loop accuracy [7], obtained in 
terms of the pole mass of the heavy quark. The pole-mass 
scheme is standard and has been used for a long time [8]. 
An alternative is to reorganize the perturbative expansion 
in terms of the running MS mass [9]. Since the correlator 
is known up to 0(a s 2 ), also the relation between pole and 
MS mass is applied to such accuracy. Figure 1 depicts the 
resulting B-meson decay constant for these two cases. 
In each case, a constant effective continuum threshold is 
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fixed by requiring maximum stability of the found decay 
constant. Thus, the constant thresholds differ for pole (so) 
and MS (so) mass schemes. Several lessons can be learnt: 

(a) In the pole-mass scheme, the perturbative series for 
the decay constant shows no sign of convergence: each of 
the LO, NLO, NNLO terms contributes with similar size. 
Consequently, the pole-mass-scheme result for the decay 
constant may significantly underestimate the exact value. 

(b) Reorganizing the perturbative series in terms of the 
MS mass of the heavy quark yields a distinct hierarchy of 
the perturbative contributions [9]. Moreover, the absolute 
value of the decay constant extracted in this scheme turns 
out to be some 40% larger than in the pole-mass case (a). 

(c) Note that, in both cases, the decay constant exhibits 
perfect stability in a wide range of the Borel parameter t. 
Thus, mere "Borel stability" is not sufficient to guarantee 
the reliability of some sum-rule extraction of bound-state 
parameters. We have pointed out this observation already 
several times [3]. Nevertheless, some authors still regard 
Borel stability as a proof of the reliability of their results. 

In the light of our above findings, we adopt in the next 
sections the OPE formulated in terms of the MS mass [9]. 



DECAY CONSTANTS OF D AND D s 

The application of our extraction procedures leads to the 
following values of the charmed-meson decay constants: 

f D = (206.2±7.3 (OPE) ±5.1 (syst) ) MeV, (4) 
f Ds = (245.3 ±15.7 (OP E)±4.5 (syst) )MeV. (5) 
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FIGURE 1. Decay constants extracted from the correlator 
given in terms of the fe-quark pole (top) and MS (bottom) mass. 



FIGURE 2. Comparison of our results for the decay constant 
fo with lattice findings. For a detailed list of references, cf. [1]. 



The OPE -related errors in (4,5) are obtained by bootstrap 
studies allowing for the variation of all QCD parameters, 
that is, quark masses, a s , and condensates, in the relevant 
ranges. We observe perfect agreement of our predictions 
with the corresponding lattice results (Fig. 2). It has to be 
emphasized that our T-dependent effective threshold is a 
crucial ingredient for a successful extraction of the decay 
constant from the sum rule (1). Obviously, the (standard) 
T-independent approximation entails a much lower value 
for the D-meson decay constant fo that resides rather far 
from both the experimental data and the lattice outcome. 

DECAY CONSTANTS OF B AND B s 

Our QCD sum-rule findings for the beauty-meson decay 
constants turn out to be extremely sensitive to the chosen 
value of the Z?-quark MS mass For instance, the 

range m b (m b ) = (4.163 ±0.016) GeV [10] entails results 
that are barely compatible with recent lattice calculations 
of these decay constants (Fig. 3). Requiring our sum-rule 
fs result to match the average of the lattice computations 
provides the rather precise value of the fo-quark MS mass 

m b (m b ) = (4.245 ±0.025) GeV. 

For this value of the b-quark mass, our sum-rule estimate 
for the B- and Z^-meson decay constants fg and fs s reads 

Sb = (193.4 ± 12.3 (OPE) ±4.3 (syst) ) MeV, 
f B , = (232.5 ±18.6 (0 PE)±2.4 (syst) )MeV. 



SUMMARY AND CONCLUSIONS 

Applying (in an attempt to improve the sum-rule method) 
our above modifications, we realize several serendipities: 

1 . The T-dependence of the effective thresholds emerges 
naturally when one attempts to render the duality relation 
exact: the dependence is evident from (1). We emphasize 
two facts: (a) In principle, such dependence on T is not in 
conflict with any properties of quantum field theories, (b) 
Our analysis of D mesons shows that it indeed improves 
decisively the quality of the related sum-rule predictions. 

2. Our study of charmed mesons clearly demonstrates 
that using Borel-parameter-dependent thresholds leads to 
lots of essential improvements: (i) The accuracy of decay 
constants predicted by sum rules is drastically improved. 

(ii) It has become possible to obtain a realistic systematic 
error and to diminish it to the level of, say, a few percent. 

(iii) Our prescription brings QCD sum-rule findings into 
perfect agreement with both lattice QCD and experiment. 

3. The beauty-meson decay constants /g,. are extremely 
sensitive to the choice of the fo-quark mass: Matching our 
QCD sum-rule fg outcome to the corresponding average 
of lattice computations provides a truly accurate estimate 
of m b {m b ), in good agreement with several lattice results 
but, interestingly, not at all overlapping with a recent very 
accurate determination [10] (for details, consult Ref. [1]). 
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FIGURE 3. Comparison of our results for the decay constant 
fs with lattice findings. For a detailed list of references, cf. [1]. 
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